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ABSTRACT 



We calculate transverse momentum distribution of direct photons from various sources 
by taking into account the initial state momentum anisotropy of quark gluon plasma (QGP) 
and late stage transverse flow effects. To evaluate the photon yield from hadronic matter we 
include the contributions from baryon-meson reactions. The total photon yield, calculated 
for various combinations of initial conditions and transition temperatures, is then compared 
with the recent measurement of photon transverse momentum distribution by the PHENIX 
collaboration. It is shown that due to the initial state anisotropy the photon yield from the 
QGP is larger by a factor of 8-10 than the isotropic case in the intermediate pt regime. It is 
also demonstrated that the presence of such an anisotropy can describe the PHENIX photon 
data better than the isotropic case in the present model. We show that the isotropization 
time thus extracted lies within the range 1.5 > r iso > 0.5 fm/c for the initial conditions used 
here. 



1 Introduction 

The primary goal of relativistic heavy ion collisions is to create a new state of matter, 
called quark gluon plasma and to study its properties through various indirect probes. Out 
of all the properties of the QGP, the most difficult problem lies in the determination of 
isotropization and thermalization time scales (ri SO and r t h erm )- Studies on elliptic flow (upto 
about pr ~ 1.5 — 2 GeV) using ideal hydrodynamics indicate that the matter produced in 
such collisions becomes isotropic with r iso ~ 0.6 fm/c [1, 2, 3]. On the other hand, using 
second order transport coefficients with conformal symmetry it is found that the isotropiza- 
tion/thermalization time has sizable uncertainties [4]. Consequently, there are uncertainties 
in the initial temperature as well. The other uncertain parameters are the transition temper- 
ature T c , the spatial profile, and the effects of flow. Thus it is very necessary to find suitable 
probes which are sensitive to these parameters. Electromagnetic probes have been proposed 
to be one of the most promising tools to characterize the initial state of the collisions [5, 6]. 
Because of the very nature of their interactions with the constituents of the system they 
tend to leave the system without much change of their energy and momentum. In fact, 
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photons (dilepton as well) can be used to determine the initial temperature, or equivalently 
the equilibration time. These are related to the final multiplicity of produced hadrons by 
isentropic expansion of the system formed in heavy ion collisions. By comparing the initial 
temperature with the transition temperature from lattice QCD, one can infer whether QGP 
is formed or not. However, it should be remembered that to characterize the initial phase 
through photons one should take into account the late stage transverse flow and early stage 
momentum anisotropy, if any. Therefore, photons could be a good probe for early stages 
of collisions, provided the observed flow effects from the late stages of the collisions can 
be understood and modeled properly. The observation of pronounced transverse flow in the 
photon transverse momentum distribution has been taken into account in model calculations 
of photon p T distribution at various beam energies [7, 8, 11, 12, 13]. 

Photons are produced at various stages of the evolution process. The initial hard scat- 
terings (Compton and annihilation) of partons lead to photon production, which we call 
hard photons. If QGP is produced initially, there are QGP-photons from thermal Compton 
plus annihilation processes. Photons are also produced in different hadronic reactions from 
hadronic matter either formed initially (no QGP scenario) or realized as a result of a phase 
transition (assumed to be first order in the present work) from QGP. In addition to that 
there is a large background of photons coming from ir° and f]° decays. The yield of excess 
photons can be obtained if this decay contribution is subtracted from the total photon yield. 
Photons from hadronic reactions and decays cannot be calculated in a model-independent 
way. The hadronic matter produced in heavy ion collisions is usually considered to be a gas 
of the low lying mesons it, p, u, rj and nucleons. Reactions between these as well as the 
decays of the p and oo were considered to be the sources of thermal photons from hadronic 
matter [5, 14, 15]. We also add the contributions from reactions involving baryons as these 
are found [16] to be comparable to that from the meson-meson reactions. 

It is to be noted that while estimating photons from QGP, it is assumed that the matter 
formed in the relativistic heavy ion collisions is in thermal equilibrium. The measurement 
of elliptic flow parameter and its theoretical explanation also support this assumption. On 
the contrary, perturbative estimation suggests the slower thermalization of QGP [17]. How- 
ever, recent hydrodynamical studies [4] have shown that due to the poor knowledge of the 
initial conditions there is a sizable amount of uncertainty in the estimate of thermalization 
or isotropization time. It is suggested that (momentum) anisotropy driven plasma instabil- 
ities may speed up the process of isotropization [18] and in that case one is allowed to use 
hydrodynamics for the evolution of the matter. However, instability-driven isotropization is 
not yet proved at RHIC and LHC energies. 

Earlier works [7, 8, 19] on photon production assume isotropy from the very begin- 
ning, i. e. Ti SO = Ti (QGP formation time). In view of the absence of a theoretical proof 
behind the rapid thermalization and the uncertainties in the hydrodynamical fits of exper- 
imental data, such an assumption may not be justified. Hence in stead of equating the 
thermalization/isotropization time to the QGP formation time, in this work, we will intro- 
duce an intermediate time scale (isotropization time, r iso ) to study the effects of early time 
momentum-space anisotropy on the total photon yield and compare it with the PHENIX 
photon data [20, 21, 22]. In the present model the space-time evolution, during the interval 
Ti < t < Ti SO , is modeled as in Ref. [23]. For the evolution from rj SO to rp (freeze-out time) 
we use (l+2)d ideal hydrodynamics. 
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Recently, it has been shown in Ref. [23] that for fixed initial conditions, the introduction 
of a pre-equilibrium momentum-space anisotropy enhances high energy dileptons by an order 
of magnitude. In case of photon transverse momentum distribution similar results have been 
reported for various evolution scenarios [24]. The model in Ref. [23] assumes two time scales: 
the QGP formation time, Tj, and the isotropization time, r iso , which is the time when the 
system becomes isotropic in momentum space. Immediately after the formation of QGP at 
Tj and Tj, the system can be assumed to be isotropic [25]. Subsequent rapid expansion of the 
matter along the beam direction causes faster cooling in the longitudinal direction than in the 
transverse direction [17]. As a result, the system becomes anisotropic with (pl 2 ) « (pt 2 ) in 
the local rest frame. At some later time when the effect of parton interaction rate overcomes 
the plasma expansion rate, the system returns to the isotropic state again (at r iso ) and 
remains isotropic for the rest of the period. 

The plan of the paper is the following. In the next section we will discuss the mechanisms 
of photon production from various possible sources and the space-time evolution of the matter 
very briefly. Section 3 is devoted to describe the results for various initial conditions and we 
summarize in section 4. 



2 Formalism 

2.1 Photon rate : Anisotropic QGP 

The lowest order processes for photon emission from QGP are the Compton (q(q) g — > q(q) 7) 
and the annihilation (qq — > g 7) processes. The rate of photon production from anisotropic 
plasma due to Compton and annihilation processes has been calculated in Ref. [26]. The 
soft contribution is calculated by evaluating the photon polarization tensor for an oblate 
momentum-space anisotropy of the system where the cut-off scale is fixed at k c ~ ^fgphard- 
Here phard is a hard-momentum scale that appears in the distribution functions. 

The differential photon production rate for 1 + 2 — > 3 + 7 processes in an anisotropic 
medium is given by [26]: 

p dN M f d 3 pi d 3 p 2 d 3 p 3 

d*xd*p ~ 2(27r)3i 2E 1 (27r) 3 2E 2 (27r)32E 3 (27r)3 /llPl ' :Phard ' Uj2lP2 ' :Phard ' U 

x (2n)H(p 1 +p 2 -p 3 -p)\M\ 2 [l±f 3 (p 3 , Phllld ,0] (1) 

where, \Ai\ 2 represents the spin averaged matrix element squared for one of those processes 
which contributes to the photon rate and Af is the degeneracy factor of the corresponding 
process. £ is a parameter controlling the strength of the anisotropy with £ > —1. f\, f 2 and 
fz are the anisotropic distribution functions of the medium partons and will be discussed 
in the following. Here it is assumed that the infrared singularities can be shielded by the 
thermal masses for the participating partons. This is a good approximation at short times 
compared to the time scale when plasma instabilities start to play an important role. 

The anisotropic distribution function can be obtained [27] by squeezing or stretching an 
arbitrary isotropic distribution function along the preferred direction in momentum space, 



m^,p h ar d ) = /n v k2 + w.iu (2) 
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where n is the direction of anisotropy. It is important to notice that £ > corresponds to 
a contraction of the distribution function in the direction of anisotropy and —1 < £ < 
corresponds to a stretching in the direction of anisotropy. In the context of relativistic heavy 
ion collisions, one can identify the direction of anisotropy with the beam axis along which the 
system expands initially. The hard momentum scale phard is directly related to the average 
momentum of the partons. In the case of an isotropic QGP, Phard can be identified with the 
plasma temperature (T). 



2.2 Photon rate : Isotropic case 

As mentioned earlier the QGP evolves hydrodynamically from r iso onwards. In such case 
the distribution functions become Fermi-Dirac or Bose-Einstein distributions. The photon 
emission rate, in isotropic case, from Compton (q(q) g — > q(q) 7) and annihilation (qq — > 
(77) processes has been calculated from the imaginary part of the photon self-energy by 
Kapusta et al. [14] in the 1-loop approximation. However, it has been shown by Auranche 
et al. [28] that the two loop contribution is of the same order as the one loop due to the 
shielding of infra-red singularities. The complete calculation upto two loop was done by 
Arnold et al. [29] and the rate is given by 

dN 1 / r „, ,„„ 1 



A(E, T) (\n[T/m q (T)} + - ]n(2E/T) + C tot (E/T)^) , (3) 



d 4 xd 3 p (2tt) 3 

where E = p and m 2 (T) = A7ra s T 2 /3 and A is the leading log coefficient given by 

m 2 (T) 

A(E, T)=2aN c J2 <& -^r 1 fn(E) (4) 
i E 

and 

C to t = C 2 ^ 2 {E/T) + C hTcms (E/T) + C^ S (E/T) (5) 

containing the dependence of the specific photon production processes. These are parame- 
terized as follows: 

- 0.3615 + 1.01 exp(-1.35£/T) 

0.548 ln[12.28 + T/E] 
(E/Tf 2 

0.133E/T 




yjl + (£/T)/16.27, 



(6) 



2.3 Photon production rate from hot hadronic matter 

First we shall consider photon emission from reactions of the type MM — > M7, where M 
generically denotes the low lying mesons. As mentioned earlier these type of reactions are 
thought to be the only sources of photons from hadronic matter and already a substantial 
amount of work has been done [5, 14, 15, 30, 31, 32] along this line. 
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Figure 1: (Color online) Medium photon spectrum, dN/d 2 pTdy, at y = for the free- 
streaming interpolating model (5 = 2) for three different values of isotropization time, Tj so , 
with initial conditions, (a) Set-I and (b) Set-Ill. 

We follow the calculations done in Ref. [32] where convenient parameterizations have 
been given for the reactions considered. These parameterizations will be used while doing 
the space-time evolution to calculate the photon yield from meson-meson reactions. The 
photon emission rate (static) from reactions of the type B M — > B7 (B denotes baryon) 
has been calculated in Ref. [16]. It is shown that this contribution is not negligible compared 
to that meson-meson reactions. To evaluate photon rate due to nucleon (and antinucleon) 
scattering from 71", p, 00, rj and a\ mesons in the thermal bath we use the phenomenological 
interactions described in Ref. [16]. 

2.4 Hard Photons 

Besides the thermal photons from QGP and hadronic matter we also calculate photons from 
initial hard scattering from the reaction of the type h,A hs — > 7X using perturbative QCD. 
We include the transverse momentum broadening in the initial state partons [33, 34]. The 
cross-section for this process can then be written in terms of elementary parton-parton cross- 
section multiplied by the partonic flux which depends on the parton distribution functions 
(PDF) for which we take CTEQ parameterization [35]. A phenomenological factor K is used 
to take into account the higher order effects. 

2.5 Space time evolution 

For any quantitative prediction of the expected total thermal photon yield the static photon 
rate, discussed in the previous section, has to be convoluted with the space-time evolution 
of the fireball. For the evolution scenario we propose the following. The system evolves 
anisotropically from to r iso where one needs to know the time dependence of Phard and 
£. We shall follow the work of Ref. [23] to evaluate the px distribution of photons from the 
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Figure 2: (Color online) Photon transverse momentum distribution with (a) Set I, (b) Set 
II and (c) Set III. The transition temperature is taken to be 192 MeV. Different lines show 
the yields for various values of Ti S0 . The data points are taken from [20, 21]. 
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first few Fermi of the plasma evolution. As the effect of transverse flow is pronounced in 
the late stages of the collisions we shall neglect this effect in the early stage. For r > r iso , 
the system is described by ideal relativistic hydrodynamics in (l+2)d [36] with longitudinal 
boost invariance [37] and cylindrical symmetry. As the system becomes isotropic at r = 
Tiso, Phard(Tiso) and T iso can be identified as the initial conditions, i. e., initial temperature 
and initial time for the hydrodynamic evolution. The time dependences of the anisotropy 
parameter £ and the hard scale Phard are taken from Ref . [23] . The initial conditions for ideal 
hydrodynamics is obtained by the conditions, 



T-ihydro / \ 

T-i = Phard(Tiso) 

hydro 



(7) 



In our calculation, we assume a first-order phase transition beginning at the time T c (phard( r c) 
T c ) and ending at r H = r d r c where r d = gqjgu is the ratio of the degrees of freedom in the 
two (QGP phase and hadronic phase) phases. We shall consider two values of the transition 
temperature T c = 192 and 170 MeV [38]. The freeze-out temperature is fixed at Tp = 120 
MeV. To cover the uncertainties in the initial conditions for RHIC energy, various combina- 
tions of Tj ,7$ consistent with the measured multiplicity (dN/dy) have been considered. We 
also vary T c to see its effect on the isotropization of the QGP. 

Therefore, the total thermal photon yield, arising from the present scenario is given by, 



dN 



d 2 p T dy 



d A xE 



dR 

d 3 p 



+ 



/ 



d A xE 



dR 

d 3 p 



(8) 



hydro 



where the first term denotes the contribution from the anisotropic QGP phase and the second 
term represents the contributions evaluated in ideal hydrodynamics scenario. 



2.6 Initial conditions and equation of state (EOS) 

To cover the uncertainties in the initial conditions for a given beam energy, we consider three 
sets of initial conditions, (I) T; = 440 MeV, n = 0.1 fm/c (II) T; = 400 MeV, n = 0.2 fm/c, 
and (III) Tj = 350 MeV, Tj = 0.25 fm/c which are consistent with dN/dy ~ 1100 measured 
at RHIC energies. The initial energy density and radial velocity profiles are taken as: 



< T ^= l + jr\ A )/S ( 9 ) 



and 



v(r h r) = v 



1 - 



(10) 



1 + e^ - ^)/ 5 . 

We also need the EOS to solve the hydrodynamic equations. Bag model type EOS has been 
used for QGP. For EOS of the hadronic matter all the resonances with mass < 2.5 GeV /c 2 
have been considered [39]. 

It is to be mentioned that in our case r iso is always less than r c and we switch on the 
transverse expansion at t 1so as the effect of transverse expansion in the very early stages is 
found to be negligible. Therefore, for r > r; so the energy density and the other thermody- 
namic variables are functions of r and r. The critical energy density corresponding to the 
quark-hadron phase transition is, thus, also a contour in the (r, r) space. 
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Figure 3: (Color online) Photon transverse momentum distributions at RHIC energies with 
initial condition Set I, for (a) T c = 192 MeV and (b) 170 MeV. 

3 Results 

We will now discuss contributions to the total photon yield due to medium photon spectrum 
from anisotropic QGP with initial conditions that might be achieved at RHIC. In what follows 
we shall consider the free- streaming interpolating model (5 = 2). The results for collisionally- 
broadened interpolating model (5 = 2/3) are described in Ref. [24]. In Fig. (1) we present the 
photon yield due to Compton and annihilation processes in the mid rapidity [9 1 = tt/2, Q 1 
being the angle between the photon momentum and the anisotropy direction) as a function 
of photon transverse momentum. Left (right) panel corresponds to Tj = 0.440 (0.350) GeV 
and Ti = 0.1 (0.25) fm/c. In estimating these results, we have used a s = 0.3. Different lines 
in Fig. 1 correspond to different isotropization times, rj S0 . We clearly observe enhancement 
of photon yield when Ti SO > Tj. The enhancement of photon yield in the transverse directions 
(y = 0) is due to the fact that momentum-space anisotropy enhances the density of plasma 
partons moving at the mid rapidity [24]. 

Next we shall consider the total photon yield from various sources. This is displayed in 
Fig. (2) for different initial conditions as described in the text. It is seen that the experimental 
data is well reproduced for all the three values of Ti SO considered here. This is because of 
the following reason. For lower values of Ti SO the initial state momentum anisotropy leads 
to lower yield as compared to higher values of r iso . But the initial temperature (phard^iso)), 
required for hydrodynamic evolution from r iso onward is higher in the former case leading to 
higher yield. These two competing effects are clearly revealed from the figure. To show that 
the presence of initial state momentum anisotropy and the importance of the contribution 
from baryon-meson reactions we plot the the total photon yield assuming hydrodynamic 
evolution from the very begining as well as with finite r iso (right panel describes the total 
contribution with and without the initial state momentum space anisotropy only for r iso = 
1 fm/c) in Fig. (3). It is clearly seen that some amount of anisotropy is needed to reproduce 
the data. We note that the value of r iso needed to describe the data also lies in the range 
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Figure 4: (Color online) Total photon yield with (a) Set II and (b) Set III with two values 
of transition temperatures, T c = 192 and 170 MeV. 

1.5 fm/c> r iso > 0.5 fm/c for both values of the transition temperatures. . 

To see the sensitivity of r iso with the initial conditions we present the results with different 
initial conditions (Sets II and III) in Fig. (4) with T c = 192 (a) and 170 MeV (b). We again 
see that the values of r iso required to fit the data are in the range of 0.5 — 1.5 fm/c and is 
almost independent of the transition temperature for a given initial conditions. 

In order to see the hydrodynamic contributions to the total photon yield we plot the 
second term of Eq.(8) in Fig. (5). The left (right) panel corresponds to Zf ydro = 318(348) 
MeV for r iso = 1(0.5) fm/c, obtained by solving 7 1 i hydro — p hard (r iso ). We also note that 
r c = 2.28(1.98) fm/c for r iso = 1(0.5) fm/c. It is found that because of the transverse kick 
the low energy photons populate the intermediate regime and consequently, the contribution 
from hadronic matter becomes comparable with that from the hadronic matter destroying 
the window where the contribution from QGP is supposed to dominate. 

4 Conclusions 

To summarize, we have calculated total single photon transverse momentum distributions 
by taking into account the effects of the pre-equilibrium momentum space anisotropy of the 
QGP and late stage transverse expansion on photons from hadronic matter with various 
initial conditions. To describe space-time evolution in the very early stage we have used 
the phenomenological model described in Ref. [23] for the time dependence of the hard 
momentum scale (phard) and plasma anisotropy parameter (£). To calculate the hard photon 
contributions we include the transverse momentum broadening in the initial hard scattering. 
The total photon yield is then compared with the PHENIX photon data. Within the ambit 
of the present model it is shown that the data can be described quite well if r iso is in the range 
of 0.5 - 1.5 fm/c for all the combinations of initial conditions and transition temperatures 
considered here. It is to be noted that the apparent hump observed in all the figures (except 
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Figure 5: (Color online) Contributions to total photon yield due to ideal hydrodynamics 
(second term of Eq.(8)) with Set II for (a) r iso = 0.5 fm/c and (b) r iso = 1 fm/c for T c = 192 
MeV. Here QM (HM) denotes quark matter (hadronic matter). 

Fig. (5)) needs to be understood and we wish to discuss it in a subsequent paper. We conclude 
by noting that the isotropization time extracted from the PHENIX photon data is within 
the limit that is required to fit the other experimental observables, such as elliptic flow at 
RHIC using ideal hydrodynamics [40]. 
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